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This paper is devoted to the method developed in Ref. |1| for 3D re- 
construction of the straight tracks in the tracking system consisting of the 
drift-chamber stereo layers. The method is based on the Hough-transform 
approach - the discrete case of more general Radon transform || - and 
takes into account both coordinates of the hit wires and drift distances not 
only for the measurements in one projection, but also in the rotated stereo 
layers. The proposed method allows one to resolve the right-left ambiguity 
and provides the accordance between vertical and horizontal projections 
of the track. 

Let the straight track of charged particle be detected by the system of 
drift chambers consisting of the cylindrical tubes placed in such a way that 
their anode wires stretched along the tube axes are parallel to the vertical 
coordinate axis Y. In this case the signals from the drift chambers give 
the information about the track projection onto the horizontal plane XZ. 
The track projection is described by x(z) = x + t x (z — z ), where x and 
t x = tan 9 X are the projection offset in z and slope, respectively. 

The single measurement from the hit wire includes its coordinates Xi, 
Zi and measured distance rj from the wire to the track. Assume that error 
of Ti measurement is uniformly distributed within the range [— c, c], where 
c » (2 -j- 3)a is a tuning parameter of the algorithm while a is a space 
resolution of the drift chamber. At each value of t x , the Hough image of 
a single measurement (xi,Zi,ri) in the space P x of projection parameters 
(xo,t x ) is given by two ranges of possible values of xq (see Fig. [l|): 
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Figure 1: The Hough image of the single measurement (xj, Zj, r^) in the space 
P x of projection parameters (xo,t x ) 



is a lower branch M d . 

For chamber operation in the proportional regime (no measurement of 
rj), the Hough image of a single hit Zi) is given in the parameter space 
P x by a range of possible values of x at each value of t x \ 
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Xi - (z { - z )t x - RJl + t 2 x ; Xi - (2^ - z )t x + RJl + tl , (3) 



where R is the tube radius. 

Assume Mi, M2, • • • , M n - are Hough images of single measurements 
from n hits belonging to the same track. In this case the Hough image of 
the track, M^ n \ is defined in the parameter space P x as an intersection of 
set of Hough images M k , k = 1, . . . , n: = QLi M k- To characterize 

the track reliability level, let us introduce the criterion J(M) which can 
be determined as number of hits having produced the Hough image M 
of the track. In particular, J(M^) = n. For the reconstruction we use 
only those M for which J(M) > J m i n , where a threshold value of J m in is a 
tuning parameter of the algorithm. If the intersection M is found, then the 
track parameters x , t x are estimated as coordinates of center of gravity of 
M in the parameter space P x . 

For 3D reconstruction, it is also necessary to determine the track pro- 
jection onto the vertical plane YZ. For this aim, the set of vertical tubes 
( "0" stereo-layers) with set of tubes rotated by the angle a around the axis 
Z ( "a" stereo-layers) can be used as it is shown in Fig. 0. The YZ track 
projection is described by y(z) = y + t y (z — z ), where y and t y = tan8 y 
are the projection offset in z and slope, respectively. 

Assume parameters xq, t x for the XZ track projection have already 
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been determined using measurements in the "0" stereo-layers and Hough- 
transform approach described above. Then, the single measurements 
(x-,Zj,rj) in "a" stereo-layers with hit-wire coordinates determined in the 
rotated system X'Y'Z (see Fig. |2|) can be used for reconstruction of track 
projection onto the plane YZ. The corresponding Hough images of single 
measurements are defined in the space P y of the parameters yo, t y by the 
following boundaries of possible values of yo at fixed t y (see Ref. for 
more detail): 

yo{ l yi x Oi l x) — 7 : tyyZi — ZQ) 

tana sin a 

Ti ± C 



± Jl + (t x cosa — tysina) 2 (4) 

sina v 

- for the drift chambers, or 

/. . \ x + t x (zi — z ) x\ 

yo{ty]Xo,tx) — " : Ty[Zi — Zo) 

tan a sin a 

R 



± J 1 + (t x cos a — t v sin a) 2 (5) 

sina v 

- for chamber operation in proportional regime. 

For the program realization of the discussed algorithm, for example, 
for XZ-projection finding, the space P x of the track parameters xq, t x can 
be treated as a discrete two-dimensional raster n x m, which is described 
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by array of its cells p(l, k) with indices I = 1, . . . ,n and k = 1, . . . , m. The 
Hough image of each hit can be consequently constructed as a strip on the 
raster: the current value of the raster cell is increased by a unit if the cell 
is located within the limits given by Eqs. (|1|), (|2]) or (|3|). After completing 
this procedure the value of each cell of the raster becomes equal to the 
number of the Hough stripes having passed the cell. Each local maximum 
on the raster exceeding some threshold J m in, can be identified with the 
XZ projection having the corresponding values of parameters Xq, t x . 

The selection algorithm of reliable tracks includes ordering the raster 
cells according to the criterion J(M) by using a special structure - the 
array of bidirectional lists determining the hierarchy of the raster cells filled 
in. The given array of the lists is a vector-column where each element is 
a pointer to the list and corresponds to a certain meaning of the criterion 
J(M) > Jmin- After filling the raster, the hierarchy of its cells is built: if 
p(l, k) > J m i n , then the point (I, k) is added into the list for J = p(l, k). 

At the first step of 3D track reconstruction, the global XZ raster and 
its hierarchy are built by using the hits only in the "0" stereo-layers. The 
further steps of the algorithm include: 

• estimating of the initial values of Xq, t x corresponding to the first 
maximum in the global XZ-raster hierarchy; 

• building the local XZ raster and its hierarchy by using the hits in 
the "0" stereo-layers within a corridor around the track projection 
defined by the initial values of xq, t x \ 

• improvement of the estimates of Xq, t x as coordinates of center of 
gravity of the local-raster maximum; 

• building the YZ raster and its hierarchy assigned to the found XZ 
projection (x ,t x ), by using hits only in the "a" stereo-layers; 

• finding of YZ projection (y Q , t y ) assigned to the found XZ projection 
(xq, t x ), by using the FZ-raster hierarchy; 

• subtraction of the hits in "0" and "a" stereo-layers, which belong 
to the 3D-track (xo,t x ;yo,t y ), and consecutive subtraction of the 
corresponding Hough strips from the global raster and its hierarchy. 

The procedure described above should be repeated iteratively until there 
remain the cells of the global XZ raster exceeding some threshold J TO j n . 

The examples of global and local XZ rasters for Monte-Carlo tracks 
are shown in Fig. 3. The small crosses on the global raster mark Monte- 
Carlo tracks on the (xq, t^-plane of the track projection parameters while 
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Figure 3: Examples of global (a) and local (b) XZ rasters 



the large crosses correspond to the values xq, t x of reconstructed pro- 
jections. Fig. 3a shows a good precision of track reconstruction by the 
Hough-transform algorithm described above. 

This algorithm has been realized in the program Htr developed for 
track finding in the PC chambers (Pattern Tracker) of the HERA-B Outer 
Tracker @]. The program Htr was integrated into the program environ- 
ment of ARTE - the general software for event processing at HERA-B - 
and tested both with Monte-Carlo and real data. The tests showed the 
stable Htr performance with average track finding efficiency of about 90% 
and rate of ghosts at the level 23% under real conditions of PC-chamber 
operation. The proposed program realization of the described algorithm 
provides the time-consuming optimization of event processing and high ef- 
ficiency of the track finding under large track-occupancy of the detector as 
well as under high level of noisy and dead channels. 

This work was done in the HERA-B software group at DESY, Hamburg. 

References 

[1] Ar. Belkov, JINR Communication P10-2001-182, Dubna, 2001. 

[2] Hough P.V.C. A Method and Means for Recognizing Complex Pat- 
terns: US Patent 3,069,654, 1962. 



5 



[3] J. Radon, Ber. Ver. Sachs. Akad. Wiss.: Leipzig, Math-Phys. KL, 
Vol. 69 (1917) p. 262. 

[4] E. Hartouni et al., HERA-B Design Report: DESY-PRC 95/01, 1995. 



6 



